N89- 24667 


CONSTRAINT ELIMINATION 
IN DYNAMICAL SYSTEMS 


R. P. Singh and P. W. Likins* 
DYNACS Engineering Company, Inc . 
Clearwater, FL 


THEME 

Large Space Structures (LSS) and other dynamical systems of current interest 
are often extremely complex assemblies of rigid and flexible bodies subjected 
to kinematical constraints. This paper presents a formulation of the govern- 
ing equations of constrained multibody systems via the application of singular 
value decomposition (SVD) . The resulting equations of motion are shown to be 
of minimum dimension. 

The motivation for this work was the development of a generic computer program 
for simulating space structures and similar electromechanical systems amenable 
to mathematical representation as a set of flexible bodies interconnected in a 
topological configuration. This representation may include closed loops of 
bodies, prescribed motion, or other constraints that may qualify as simple 
monholonomic . The equations of motion appropriate for a set of flexible bodies 
in an open loop configuration appear in Refs. 1, 2. A computer program 
(TREETOPS) developed to simulate the dynamic response of flexible structures 
in a topological tree configuration is described in Ref. 3. The SVD technique 
of the present paper is being incorporated in an extension of the TREETOPS 
program that permits application to constrained systems. This extension 
permits direct use of the dynamical equations for the less constrained system 
in Refs. 1, 2, with augmentation by kinematical constraint equations and re- 
duction of dimension by SVD. 

Basically, there are two conceptual approaches to solving the equations of 
motion of such systems. (1) One can introduce unknown forces and torques at 
the interfaces between constrained bodies (often accomplishing this symbolical- 
ly with Lagrange multipliers) , and then solve the dynamical equations simul- 
taneously with the constraint equations to determine the constraint forces and 
torques as well as the kinematical variables. Ref. 4. (2) Alternatively, one 

can use the constraint equations to reduce the dimension of the system of 
dynamical equations to be solved by partitioning generalized coordinates. Refs. 
5, 6. Techniques presented in Refs. 4, 5, 6 may encounter numerical singular- 
ities. Also, systems undergoing large motion may present problems of 
inconsistency in the constraints such as three-dimensional loops during the 
system motion becoming two dimensional or one-dimensional loops. In what 
follows, the SVD method will be shown to avoid mathematical singularities. 
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Singular Value Decomposition: Orthogonal decomposition of an mxn matrix L by 
singular value decomposition is closely related to the eigenvalue-eigenvector 
decomposition of the symmetric positive semidefinite matrices L L and LL . 

Let r<m be the rank of L. Then there are orthogonal matrices U and V of order 
mxm and nxn rexpectively such that 


T 

U LV = 


Z 0 
0 0 


( 1 ) 


= s 


wher e E — drag (A , Af.**.,A ) and A^A^...^A ^ 0 • 

1 2 r 12 r 

The diagonal elements of the decomposition are called the singular values of 
the matrix L. The singular values are unique, although U and V are not. 


It is easy to verify that 

V T L T LV = diag (E 2 , 0) (2) 

2 9 T 

Thus (A*,....,A ) must be the nonzero eigenvalues of L L arranged in the 

descending order and the requirement that A. be nonnegative completely 

T T 

determines the A^. The eigenvectors of L L are the columns of V. If L L 

has a multiple eigenvalue A 2 >0, the corresponding columns of V may be chosen 
as an orthonormal basis for the space spanned by the eigenvectors corresponding 
to A 2 . 


From eq. (1) 

T 

L = USV 


(3) 


Now with proper partitioning of U and V eq. (3) can be expressed as 



(4) 


From the above one obtains 
-1 

Uj, = LV X Z (5) 

Thus once V is chosen U is obtained by eq. (5). The matrices U and V may 
II ^ T ^ 

be any matrices with orthonormal columns spanning the null spaces of L and L, 
respectively. It is worthwhile to mention that the null space of L is the 
space of all vectors x such that 
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Lx = 0 


( 6 ) 


With the orthogonal decomposition given by eq. (3), an nxm matrix L , called 
the pseudoinverse of L r is defined by 



L + is uniquely defined by L; it does not depend on the particular orthogonal 
decomposition of L. 

Application of SVD to Dynamical System with Constraints : Let q = q , . . . . ,q n 

comprise a set of generalized coordinates that fully defines the configuration 
of the dynamical system. The equations of motion of the system can be written 
as 


Mq = F(q,q,t) (8a) 

where the elements of nxn matrix M are functions of q's and the inertia pro- 
perties of the system; the elements of nxl column vector F are functions of q's, 
their time derivatives q's and applied forces (moments) on the systems. If 
the generalized coordinates are related by constraint equations then they are 
not independent and the right hand side of eq. (8a) will also include the non- 
working forces of constraints. Let the unknown constraint forces be denoted 

Q 

F . Now for the general case of constrained dynamical system, eq. (8a) takes 
the following form 

Mq * F + F c (8) 

Suppose however that the constraint equations can be written as 

Aq = B (9) 


where A is of dimension mxn (m<n) and B is an mxl column vector. 

Holonomic constraint equations can always be placed in the form of eq. (9) and 
nonholonomic constraints in the class called Pfaf f ian or simple have this 
structure also. 

If the rank of matrix A is r<m then r of the kinematical variables in q are 
related by eq. (9) and there are only n-r independent generalized coordinates. 
In other words the dynamical system possesses n-r degrees of freedom. 

The SVD of the mxn matrix A provides 

A = USV T (10) 

The orthogonal matrices U and V (of dimension mxm and nxn, respectively) are 
partitioned as 
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where and are respectively mxr and nxr matrices; and are respec- 
tively mx(m-r) and nx(n-r) matrices. Note that r is the rank of A. 

Because AV 2 = 0, eq. (9) is satisfied by 

q = A + B + V z (13) 

for any vector z, A + is the pseudoinverse of A. We shall refer to z as the 
reduced set of (n-r) coordinates. 


Differentiation of eq. (9) with respect to time yields 

M • • • 

Aq = -Aq+B 

(14) 

(13) express q in terms of z as 

+ V 2 z (15) 

(13) or eq. (15) that maps the n kinematic variables q (or q) 

to n-r variables z (or z) . Thus a consistent set of equations of motion in z 
is given as 

T •• T T c T + 

V 2 M V = V 2 F + V 2 F - V 2 MA B' (16) 

The coefficient of z is a symmetric positive definite matrix with the charac- 
teristic of an "inertia matrix" for the reduced set of coordinates z. 

With the Lagrange multiplier method, F° is established via (see Ref. 4) 

F C = A T 0t (17) 


or, Aq = B 1 

Following eq. 

.. + , 
q = A B * 

Note from eq. 


where a is the column vector of Lagrange multipliers. 

T 

Premultiply eq. (17) by V 2 to obtain the following 

T c T T 
V 2 F = V 2 A a 


- (AV 2 ) a 


= 0 


(18) 


540 



Thus it is seen that the nonworking constraint forces make no contribution to 
the equations of motion (eq. (16)) and need not be recorded. 

Employing the transformations given by eqs. (13) and (15) , the minimum 
dimension governing differential equations of motion are given by 

T - T T + 

V MV z = V F - V MA B 1 (19) 

2 2 2 2 


and 


4 - A + B 4 - V 2 z (20) 

This method eliminates the forces of constraints which when included serve not 
only to enlarge the dimension of the dynamical system but also quite often 
introduce computational problems. 
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